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ABSTRACT 



The aftermath of supernova (SN) 1987A continues to provide spectacular 
insights into the interaction between a SN blastwave and its circumstellar en- 
vironment. We here present 36 GHz observations from the Austraha Telescope 
Compact Array of the radio remnant of SN 1987A. These new images, taken in 
2008 Apr and 2008 Oct, substantially extend the frequency range of an ongo- 
ing monitoring and imaging program conducted between 1.4 and 20 GHz. Our 
36.2 GHz images have a diffraction-limited angular resolution of 0'/3 — 0'.'4, which 
covers the gap between high resolution, low dynamic range VLBI images of the 
remnant and low resolution, high dynamic range images at frequencies between 1 
and 20 GHz. The radio morphology of the remnant at 36 GHz is an elliptical ring 
with enhanced emission on the eastern and western sides, similar to that seen 
previously at lower frequencies. Model fits to the data in the Fourier domain 
show that the emitting region is consistent with a thick inclined torus of mean 
radius 0'/85, and a 2008 Oct flux density of 27 ±6 mJy at 36.2 GHz. The spectral 
index for the remnant at this epoch, determined between 1.4 GHz and 36.2 GHz, 
is q; = —0.83. There is tentative evidence for an unresolved central source with 
flatter spectral index. 

Subject headings: supernovae: general, supernovae individual: (SN1987A), supernova 
remnants, supernova remnants individual: (SN1987A), radio continuum: general 
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Introduction 



The type IIP supernova 1987A is the most well-observed and published supernova 
in existence. As the brightest supernova seen since the invention of the telescope, it 
has been observed right across the electromagnetic spectrum from radio through to 
gamma rays. It was also the first supernova for which we have observations of the 



progenitor, as well as th e first and only supernova thus 



a neutrino burst event OHirata et al.l |198yI : iBionta et al. 



'ar to have been associated with 



19871 ). Neutrino observations 



from SN1987A have bee n used to pu t a firm foundation on theories of core collapse and 



neutron star formation ( Jankj 



19971 ) as well as giving a time of core collapse at 1987 



February 23, 07:35 UT. Shortly after core collapse an expanding shockwave tore through 
the progenitor's stellar envelope. The resulting sho ck breakout resulted in a UV flash that 



ionized material in the circumstellar environment ( Luo fc McCray 



199ll ). illuminating a 



central equatorial ri ng and two outer ring s that we now know are part of a much larger 



hourglass structure (ISugerman et al 



20051 ). As yet, the mechanism for ring formation is 



not properly understood, however two competing theories, the interacting wi nds model 



( iBlondin fc Lundqvist 



binary merger model (jPodsiadlowski et al. 



1993 



Martin fc Arnett 



1995 



Tanaka fc Washimi 



2002) and the 



20071 ) offer plausible mechanisms for shaping 



the rings. At present the shock is transiting the equatorial ring and continues to provide 
insights into the hydrodynamics of a unique young supernova remnant as followed by X-ray, 
optical and radio telescopes. 

In the radio SN 1 987A was imaged as soon as a few days after detection. On 1987 
Feb 25, 



Turtle et al. 



(119871 ) detected emission at 843 MHz with the Molonglo Observatory 
Synthesis Telescope (MOST) along with detections at 1.4, 2.3 and 8.4 GHz by an 
interferometer formed between the Parkes 64 m and Tidbinbilla 34 m antennas. Emission 



at 843 MHz peaked around 130 mJy on 1987 Feb 26-27. Just over five days after the 
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neutrino burst (1987 Feb 28), IJauncey et al.l (Il988l ) used VLBI techniques to resolve the 
expanding radiosphere. Continued observations saw the su pernova fade frorn . its intial radio 



19871 ). From this 



burst to almost undetectable levels within a few months (iTurtle et al. 
point on, radio emission from SN1987A deviated radically from the trend of typical type 
IIP supernova light curves. Radio emission fr om the remnant returned in 1990 with a 



re-detection at 843 MHz by MOST on July 5 flBall et al. 



200ll ) and the six 22m antenna s 



of th e Australia Telescope Compact Array (ATCA) a month later (jStaveley-Smith et al. 



19921 ). This has been interpreted as a collision between the expand ing shock front and an 



ioniz ed hydrogen region about half the distance to t 
1995 ) . Since then, an ongoing observing campaign ( 



2002 



le optical ring 



2005 



Staveley-Smith et al. 



2007 



Gaensler et al. 



Gaensler et al. 



Chevalier fc Dwarkadas 



1997 



Manchester et al. 



20071 ) has followed radio emission 



from the developing remnant at frequencies ranging from 1.4 to 20 GHz. This saw the 
radio emissi on rise at monotonic ally increasing rate as the shockwave began its crossing 



of the ring (jZanardo et al 



20091 ). Such behaviour puts the remnant of SN1987A in the 
somewhat rare category of young and brightening supernova remnants, along w ith the 



Galactic supernova remnant Gl-9+0.3 (IReynolds et al. 



2008 



Murphy et al 



20081). 



High precision images of the remnant at high resolution are important for studying 
the morphology and e volution of the b lastwave interaction. Previous estimates of the 



optical ring geometry (jPlait et al 



19951 ) have shown that it is an ellipse with semi-major 



axis 0'.'858 ± O'.'Oll and semi-minor axis 0'.'621 ± O'.'Oll, thus requiring very high resolution 



imaging. Cameras ab oard the Hubble Space Te 



at O'.'Oa - 0'.'2 pixel"^ fiMcMaster fc Birettal 120081 1 which has enabled detailed imagmg 



escope (HST) are capable of imaging 



of the structure of t he equatorial ring as well as the development of the reverse shock 



(iMichael et al 



20031 ). Pixels on the ACIS detector aboard the Chandra X-ray Observatory 



register half the total encircled energy of incident photons within a circle of diameter 0'.'8 
[half power beam diameter (HPD)], although this can be improved to around 0'.'6 HPD with 



subsampling techniques (iTsunemi et al. 



200 ID . In 2007 October, 



Tingay et al 



(120091 ) used 



an e-VLBI array formed between the ATCA, Parkes, and Mopra (22m) antennas to image 
the remnant of SN1987A at 1.4 GHz with a restoring beam of 85x168 milh-arcseconds. 
However the hmited range of baselines made such an array insensitive to details larger than 
0"A. In the radio the ATCA is capable of resolving details as small as 0'.'2 at around 40 
GHz, which nicely fills in the gap left by the VLBI observation and images at the lower 
frequencies. Imaging at higher frequencies around 100 GHz is also possible, although this 
capability is enable to for only five of the six antennas, cutting the number of baselines 
by a third as well as reducing the maximum baseline length to 3 km. This translates to 
roughly the same resolution (approximately 0'.'2) as the 40 GHz band. Unfortunately at 
higher frequencies atmospheric opacity degrades t he image with approximately an order 



of magnitude higher effective system temperature (IBrooks 



20071 ). With an aim to obtain 



the highest res olution iraages c urrently possible with the ATCA, the ongoing imaging 



observations of 



Ng et al 



(120081 ) were extended to include frequencies around 36 GHz. In 
this paper we present the latest 36 GHz data from the recent extension to these ongoing 
observations. 



2. Observations & reduction 

Our observations of the supernova remnant were made with the ATCA|^ at center 
frequencies of 34.88 and 37.44 GHz on 2008 April 25 and 2008 October 7,8,12 using antenna 
configuration 6A with multifrequency imaging over 2 x 128 MHz bandwidth in all Stokes 



^The Australia Telescope Compact Array is part of the Australia Telescope which is 
funded by the Commonwealth of Australia for operation as a National Facility managed by 
CSIRO. 
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parameters. Atmospheric phase stabihty at the ATCA was monitored with two 1.8 m 
antennas spaced 230 m apart designed to measure the path noise of a 30 GHz signal from 
the geostationary Optus-B3 communications satelhte. A path noise of below 300 /im is 
considered good for millimetre-wave astronomy at the ATCA. Table [1] shows observing 
conditions and un-fiagged uv coverage for each observation. 

Table 1: Observing conditions for 36.2 GHz observations of SN1987A. 



Date 




Path noise 
(/im) 


Integration time 
(hr) 


2008 


Apr 25 


334 ± 133 


11.3 


2008 


Oct 7 


131 ±22 


7.4 


2008 


Oct 8 


88 ±10 


3.9 


2008 


Oct 12 


184 ± 60 


9.8 



The 2008 Apr 25 observation was hampered by poor weather, consequently around 3 
hours of uv-data was discarded. The secondary flux calibrator PKS 0637-752 was observed 
for 2 minutes in hourly intervals and the phase calibrator PKS 0530-727 was observed for 
1.5 minutes per 6 minutes integration time on source. A flux density of 3.93 ± 0.06 and 
3.83 ± 0.06 Jy at 34.88 and 37.44 GHz was obtained for the 2008 October observations 
of PKS 0637-752 by bootstrapping an observation of the primary flux calibrator Uranus 
at similar elevation 35°) from the night of 2008 October 8. The flux density of PKS 
0637-752 for the 2008 April 25 observation was similarly determined to be 3.71 ± 0.07 and 
3.61 ± 0.08 Jy using a later 2008 May 25 observation of both the calibrator and Uranus 
(at the shortest baseline). The phase center was located about 5'.'5 south of the remnant to 
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avoid defects (ITaylor et al.lll999l ). We used the software package Miriad (ISault et al.lll995l ) 
to perform image reduc tion and analysis. Continuum images were made with robust=0.5 



weighting (iBriggs 



19951 ) as this represented the best compromise between resolution and 
signal-to-noise ratio. Natural weighting (robust=2.0) was used in both the flux and 
polarisatio n measurements in order to preserve as much signal as possible. The CLEAN 



algorithm (iHogbom 



1974 



Clark 



1980l ) was employed to deconvolve the resulting continuum 
images. Table [2] gives the size and position angle of the restoring beam, as well as the rms 
noise of the resulting image in each case. 



Table 2: Image parameters for the 2008 observations at 36.2 GHz. 



Date 


Robust 


Restoring beam 


Position angle 


Rms noise 






(") 


n 


(mJy beam~^) 


2008 Apr 25 


0.5 


0.34 X 0.24 


45.0 


0.10 




2.0 


0.38 X 0.28 


44.6 


0.10 


2008 Oct 7-12 


0.5 


0.37 X 0.22 


-7.7 


0.06 




2.0 


0.43 X 0.25 


-8.3 


0.06 


2008 Apr-Oct 


0.5 


0.33 X 0.24 


-1.3 


0.06 




2.0 


0.38 X 0.27 


-1.0 


0.06 



Shown in Figure [T] are the deconvolved Stokes-I continuum images of the 36.2 GHz 
observations. Figure [H^a) is the image from the 2008 April 25 observation, Figure [11(b) is 
the combined 2008 October image, and Figure [T](c) is the combined 2008 April and 2008 
October images. Note that the approximately north-south extensions prominent in the 
images are probably artifacts caused by insufficient uv-coverage and increased phase errors 
when observing at low elevation. 
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Fig. 1. — Continuum images of SN1987A made with the ATCA at an average frequency of 
36.2 GHz and robust=0.5 weighting: (a) 2008 April 25; (b) 2008 combined October data; 
and (c) 2008 combined observations. The 2008 April, October and combined images have 
an rms noise of 0.10, 0.06 and 0.06 mJy beam~^ respectively. The linear greyscale on each 
image ranges from to 3 mJy beam~^. The black ellipse at lower left is the beam as given 
in Table E) 
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For purposes of comparison, we also obtained a contemporaneous observation at 
8.6 GHz on 2008 Oct 9-10 as well as our regular rn o nitori ng; at 1.4-8.6 GHz. The 8.6 GHz 



observation was reduced as described in 



Ng et al. 



fl2008h . and the 1.4-8.6 GH z flux 



monitoring observations and reduction are as described in (jZanardo et al. 



Manchester et al 



(120021 )]. 



20091 ) [see also 



3. Results and discussion 

3.1. Flux density and spectral index of the remnant. 

A flux density estimate was obtained for the naturally weighted (robust=2.0) 2008 
April and October images. The amplitudes and phases of the April observation were 
affected by prevailing weather conditions, making the 2008 April 25 flux density estimate 
unreliable. A flux density of 27 ± 6 mJy was measured for the combined October image by 
imaging using natural weighting and taking a cut 6a above an rms noise of 60 /iJy beam"^. 
The flux density of the 1.4, 2.4, 4.8, and 8.6 GHz 2008 October monitoring observations 
was measured to be 414 ± 12, 280 ± 8, 173 ± 6, and 112 ± 6 mJy, respectively. This yields 
a spectral index a = —0.71 ± 0.01. When used to predict the 2008 October flux density 
at 36.2 GHz, the fltted power spectrum predicts a flux density of 41 ± 6 mJy. The 2008 
October flux density of 27 ± 6 mJy therefore lies 1.3cr lower than the power law flt to the 
other frequencies. If we include the 2008 Oct observation in the flt, we obtain a spectral 
index of a = —0.83, where oc ly". A comparison of the flts is shown in Figure [2l 



3.2. Modeling in the Fourier domain and a limit on a central source. 



We modeled the 2008 Oct 36.2 and 8.6 GHz d ata in the Fourier domain using a 
truncated shell model described by iNg et al.l (120081 ). The shell has a radius R and a 
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Fig. 2. — Spectral index of SN1987A as determined by a power-law fit to the flux densities of 
the 2008 October monitoring observations both with and without the 36.2 GHz observation. 
The spectral indices from the fit are a = —0.71 for the solid line and a = —0.83 for the 
dashed line. 



half-opening angle 9 (i.e., 9 = 90° and 0° corresponds to a spherical shell and an equatorial 
ring respectively), and it is inclined at 43.4° south of the line of sight at a position an gle 
-7.6°(or 7.6°clockwise) from north, as determined from optical observations (jPunI 120071 ). We 
modeled the east-west asymmetry by a line a r grad ient in the equatorial plane, more details 



of the modeling can be refered to 



Ng et al. 



(120081 ). The best fit model parameters at the 
two frequencies are shown in Table [3] and displayed graphically in the top row of Figure 
[3l Both models are roughly consistent, each with 30% greater flux in the eastern lobe. 
The 8.6 and 36.2 GHz fitted shell widths do not agree at the 3a level. If this is a physical 
effect, it may be due to a positional dependency of the distribution of particle energies in 
the forward shock. The small shell width characteristic of the fits to low frequency radio 
emission may correspond to a proportionately larger population of newly swept up particles 



at the forward shock. Ho wever she , 
by the fit, as discussed in iNg et al. 



1 thick ness 5 is a parameter that is not well constrained 



(120081 ) and demonstrated in Table 2 of that paper. 



Table 3: Best fit truncated-shell models of the 2008 Oct 36.2 GHz and 2008 October 8.6 GHz 
observations. 



Frequency Flux density Mean shell radius Shell width Half-opening 

f R 5 angle 9 

(GHz) (mJy) f) (%) (°) 

36.2 GHz 32.7 0.850± 0.060 16 ±4 29.2 ±2.0 

8.6 GHz 108.0 0.892 ±0.002 0.04 ±2.6 35.6 ±0.6 
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At 36.2 GHz, the model was subtracted from the data to produce a residual image as 
shown in Figure |l](c), alongside the imaged data and imaged model in Figures Hl^a) andHI^b). 
Features present or absent in the fitted model appear as negative or positive components 
in the residual image. Within the remnant in the residual image, the highest positive peak 
is +0.2 ± 0.1 mJy beam~^ and the lowest negative peak is —0.2 ± 0.1 mJy beam~^ with 
an average difference of zero. The amplitude of the peak, the maximum difference between 
the real and model image, represents an upper limit on the flux density from a possible 
central source. A more quantitative measurement was made by adding an additional 
point source to the truncated torus model and fltting the model to the 2008 October 
uv-data at 36.2 GHz. The resulting flux density of the fltted point source in the model is 
0.3 ± 0.2 mJy, which is consistent with the above est i mate . It is also consistent with the 



0.3 mJy upper limit at 9 GHz, estimated in 



Ng et al. 



(120081 ). Possible reasons for this extra 



point source flux include: an energized population of accelerated electrons due to projected 
shell emission at high latitudes away from the equatorial ring, a source within the central 
region of the remnant, noise, or possibly phase errors in the 36.2 GHz image. Assuming 
0.3 mJy as the flux density of a possible central source, we e stimate a 36 G Hz spectral 



luminosity of 1 x 10^^ Watts Hz ^ at a distance of 51.4 kpc (jPanagia 



20031 ). This value 



is signiflcantly under- 
supernova remnants. 



uminous when co mpared to some central sources of other known 



Page et al 



( 120071 ) obtain a flux density of 322 Jy at 33 GHz and 
299 Jy at 40 GHz from WMAP observations of the Crab Nebula. Assuming an interpolated 
flux density of 310 Jy f or 36.2 GHz g ives a luminosity for the Crab of 1.5 x 10^^ W Hz~^ 



at a dista nce of 2 kpc ( 



Trimble 



SN1986J (IBietenholz fc Bartel 



19731) . Observations at 22 GHz of the young remnant of 



20081 ) show a flux density of 4 mJy in a central region of flat 



spectral index. Extrapolated to 36.2 GHz, this corresponds to a flux density of 3 mJy or a 
luminosity of 4 x 10^^ W Hz~^ at a distance of 10 Mpc. 
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3.3. Spectral index image 



SN1987A has been 'c 



around 14.7 solar masses (INomoto et al 



assified' as a Type HP supernova with a progenitor mass of 



19891). We know from neutrino observations that 



a large number of neutrons were formed in the core collapse. However it is still uncertain 
if a neutron star su rvived the explosi on. Searches for an associated pulsar have been 



unsuccessful so far (Manchester 



20071), although they only rule out the existence of a 
relatively high-luminosity pulsar. Another possible way to find evidence of a surviving 
neutron star is to look for a pulsar wind nebula. Pulsar wind nebulae (PWNe), are 
clouds of non-thermally radiating particles that are excited by an enclosed compact 
central source. PWN e have been detected in about ten percent of all known core collapse 



supernova remnants (IKaspi fc Helfand 



2002), with the Crab nebula and 3C58 being some 



of the the best-known examples. In the radio such nebulae are characterised by polarised 
sy nchrotron emission and typically have a flat spectral index of ranging between -0.3 and 



(iGaensler fc Slane 



20061 ) . Detection of fiat spe ctral index emission n ear the explosion 



2004 ) ■ would be a good 



site of SN1987A, as was the case with SN1986J (iBietenholz et al. 
indication that a neutron star has survived the explosion. A search for a possible central 
fiat spectrum source was undertaken by generating a spectral index image using the 2008 
Oct data at 36.2 GHz and the 2008 October data at 8.6 GHz [discussed in Section |2]. As 
the synthesised beam for the 8.6 GHz data is 0'.'4, the 36.2 GHz data were re-imaged with 
a 0'.'4 restoring beam using the MAXEN algorithm instead of CLEAN for consistency with 
the 8.6 GHz image. Both the 8.6 GHz and the re-imaged 36.2 GHz images are shown in the 
middle row of Figure [31 As a precaution against possible false spectral index information 
generated by the imaging process, we also produced a comparison spectral index image 
made from the Fourier models in the top row of Figure [3l The models were inserted into the 
visibility data with a suitable position offset to locate it at the same angular distance from 
the phase center but on the opposite side of the phase center from the remnant. Following 
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imaging and deconvolution, the resulting model images, shown in the bottom row of Figure 
131 have approximately the same noise, flux density, and restoring beam characteristics as 
their real counterparts. In Figure E] we present the two spectral index images. Figure 
[5]^a) was formed by comparing the 2008 Oct 8.6 GHz and re-imaged 2008 Oct 36.2 GHz 
observations, and Figure [5](b) was formed by comparing the 2008 Oct 8.6 GHz and 2008 
Oct 36.2 GHz model images. The real spectral index image has been masked using the 
same mask used to determine the flux of the 2008 October CLEAN image at 36.2 GHz. A 
similar mask has also been applied to the model spectral index image. 

Both spectral index images share the same broad characteristics, regions of spectral 
index, a ~ —0.8 at the two lobes and regions of steep spectral index, a ^ —1.5 in the 
interior. There are two notable flat spectral index features in the real image. The feature 
at the northern rim is caused by the spurious north-south extensions mentioned in Section 
[21 The feature around 0'.'2 east of the progenitor position, which has a maximum spectral 
index of a ~ —0.7 ± 0.3 and encompasses an area of approximately one beam (width 0.4"). 
The position of the flattest spectral index in this feature is 05^35'"27'.991, -69°16'11'.'13 
(J2000) with around 50 mas positional uncertainty. However a smaller flat spectral index 
feature is also present in the model counterpart. This feature is largely due to the imaging 
process rather than noise or differences between the models. Follow-up observations will 
be needed to determine if the feature in the real spectral index images corresponds to an 
actual source of emission in the central region emission not fitted by our Fourier models. 



3.4. Polarisation and magnetic field direction. 



Magnetic field direction in supernova remnants can be divided into three categories: 
tan gled, radial and t angen tial. In the atlas of supernova remnants compiled by iMilnd (119871 ) 
and iFiirst fc ReichI (120041 ). 27 out of the 70 remnants surveyed showed polarisation. The 
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absence of polarisation in the presence of synchrotron emission indicates a tangled magnetic 
field and/or extreme rotation measures along the line of sight. Of those remnants which 
did exhibit polarisation, such as Cassiopeia A and Tycho's remnant, most under 10,000 
years old had radial magnetic field directions across the shock front. Ray leigh- Taylor (R-T) 
instabilities have been implicated as the mechanism for the radial stretching of the magnetic 
field in young supernova remnants, however magnetohydrodynamical simulations have 



yet t o produce R-T instabilities that reach across the forward shock (IBlondin fc Ellison 



2OOII ). It is anticipated that higher shock compression ratios generated by the influence of 



shock-acce^ 



erated partic 



the shock (ISchure et al. 



es w ill enable the generation of radial magnetic fields that cross 



20081 ). Older supernova remnants such as Vela exhibit a tangential 
magnetic field, interpreted as the accumulation of an interstellar magnetic field at the 
ageing blast front. 

A search was made for linear and circularly polarised emission in the remnant 
of SN1987A by constructing polarisation intensity images with natural weighting. No 
persistent linear polarisation in Stokes-Q was detected above a 3a threshold of 0.6 and 
0.24 mJy beam"^ for the 2008 April and 2008 October observations. A similar technique 
was employed to search for circular polarisation in Stokes-V, however the same null result 
was obtained for a 3a threshold of 0.6 and 0.24 mJy beam~^. Possible causes for the 
de-polarisation of the radio emission include a tangled magnetic field in the emitting region 
and/or extreme Faraday rotation measures on the order of ~ 50000 rad m~^. 



3.5. Radial profile and astrometry 



The optical position of the supernova at day 1278, as measured by a frame-tie to 
the Hipparcos and VLBl reference frames is 05'^35'"27^9 68, -69°16ll^^09 ('J20 00), with 



an error of 30 milli-arcseconds (mas) in each coordinate ([Reynolds et al. 



I995I). At day 
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7899, corresponding to 2008 Oct 9, this position will have mo ved due to the proper mo tion 



Kallivayalil et al 



fl2006h is 



of the remnant. The proper motion of the LMC according to 
+2.03 ± 0.08 mas yr~^ east and +0.44 ± 0.05 mas yr~^ north. The remnant also experiences 
a 0.3 ± 0.02 mas y r~^ motiori north due to its estimated 65 + 5 km s^^ bulk rotation speed 



around the LMC (iKim et al. 



19981 ) at a distance of approximately 5 kpc from the LMC's 
kinematic center (assuming a distance of 50kpc to SN1987A). If corrected for both effects, 
then the progenitor should, at day 7899, be situated 37 + 1.5 mas east and 13 + 1.3 mas 
north of Reynold's position, assuming an LMC distance of 50 kpc. This predicted position 
corresponds to the coordinates of 05'^35™27'^.975, — 69°16'11'.'08 (J2000) with an uncertainty 
of 30 mas in each coordinate. We tested this prediction by fitting an ellipse to points that 
make up a ridgeline of the radio ring in the combined 2008 April-Oct 36.2 GHz image. 
Points along the ridgeline in the radio image were located by propagating rays away from 
the original progenitor position. The point of maximum flux density along each propagated 
ray was determined to be part of the ridgeline. We fitted an ellipse to both fiux-weighted 
and non fiux-weighted ridge line points and tabulated the fitted ellipsoids in Table HI We 
see that fitting an ellipse to the 36.2 GHz October data in the image plane produces a 
semi-major axis that is slightly smaller but consistent with the mean radius obtained by 
fitting a truncated-shell model to visibility data in the uv plane. 



In Figure [6] the fiux-weighted ellipse and ridgeline points are overlaid atop the 
combined 2008 Apr-Oct 36.2 GHz image along with the two fitted ellipse centroids 
and the original progenitor position. The fiux-weighted ellipse centroid was located 
at 05^35'^27^987, -69°16'11'.'115 (J2000) and the un-weighted centroid was located at 
05^35"^27^990,-69°16'll'.'110 f J2000). Becaus e we h ave used the VLBI position of the 



phase calibrator 0530-727 from .Reynolds et al. 



(119951 ) ■ the estimated positional uncertainty 



of the fitted centroids above is about 30 mas in each coordinate with an additional ~ 15 — 17 
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Table 4: Parameters of the flux-weighted and flux un-weighted elliptical fits. 



Weighting 


RA offset 1 


Dec offset^ 


Semi-major axis 


Semi-minor axis 


Position angle^ 




(mas) 


(mas) 


(") 


(") 


n 


Weighted 


+72 ± 48 


-37 + 46 


0.79 


0.68 


74.7 


None 


+79 ± 48 


-32 + 46 


0.80 


0.66 


70.6 



Notes. ^Offsets are the angular distances of the fitted ellipse centroids from the 
predicted position of 05^35"^27^975, -69°16'11'.'08 (J2000) 
^Position angle is defined North through East. 

mas uncertainty in the fit. As shown in Figure O the offset of the flux-weighted centroid is 
72 + 48 mas east and 37 + 45 mas south of the predicted position, and the flux-unweighted 
centroid is 79 + 48 mas east and 32 + 46 mas south of the predicted position. Given that 
the positional uncertainty in the predicted position is about 30 mas, we conclude that 
the two fitted centroids and the predicted position all lie within la of each other. If a 
hydrodynamical effect following the explosion is the cause of the centroid offset from the 
predicted position, we estimate that the (day 7899) upper limit for any departure from 
symmetry for the eastern lobe, estimated as approximately twice the offset plus uncertainty, 
is at most 174 mas east and 64 mas south. The position of the possible flat spectral index 
feature in Figure [S](a) was offset 124 mas east and 32 mas south of the predicted position, 
which places it within 60 mas of the fitted centroids. 
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Fig. 3. — Fourier models and images from the 2008 October observation at 8.6 GHz and 
36.2 GHz. The left column contains 8.6 GHz images, the right column 36.2 GHz images 
obtained using the same MAXEN deconvolution algorithm. In the top row are the truncated- 
shell models, the middle row contains the real images. The bottom row contains model 
images formed by inserting the model visibility data, generated from the top row, into the 
observed visibility data from the middle row and running them through the same imaging 
pipeline. The model images (bottom row) thus include noise and have the same uv-coverage 
as the real observations (middle row). All images have been restored using an identical 0'.'4 
beam and have a linear greyscale that ranges from 0-13 mJy for the 8.6 GHz images and 
0-3.5 mJy for the 36.2 GHz images respectively. 
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Fig. 4. — (a) 2008 October 36.2 GHz image of SN1987A; compared with (b), the model from 



the Fourier-domain fitting flNg et al.ll2008l ): and (c) the imaged residual image between the 



image and model. The real and model mages are displayed on a linear greyscale ranging 
from to 4 mJy beam~^, and the residual image has a greyscale ranging ±0.3 mJy beam~^. 
The rms noise of the residual image is estimated to be around 90 /iJy. 
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A RA (arcsec) A RA (arcsec) 

Fig. 5. — Spectral index images formed by comparing the super-resolved 2008 Oct 8.6 GHz 
and the re-imaged 2008 Oct 36.2 GHz images of SN1987A with the same MAXEN decon- 
volution algorithm and 0'.'4 restoring beam, (a) is the real spectral index image and (b) is 
the model spectral index image generated from the imaged Fourier models. The progenitor 
position is shown white cross in each image. 
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Radial profiles of the remnant were made by propagating rays away from the averaged 
position of the ellipse centroid [obtained from the ridge-finding algorithm] and plotting the 
resulting profiles, shown in Figure [71 In particular, note that the flux density within the 
projected radio ring on the image plane is always greater than 6a above the rms noise of 
0.06 mJy beam~^, which indicates that emission is being produced along the line of sight, 
either from high latitude regions of the expanding shell, or from within the central region 
of the nebula itself. 



3.6. Comparison with observations at other wavelengths 



In order to see how the 2008 Oct 8.6 GHz and 2008 Oct 36.2 GHz images compared 
with observations at other wavelengths we obtained optical Hubble Space Telescope (HST) 
and X-ray Chandra images from their respective public archives. The HST image was taken 
on 2006 Dec 6 using filter F625W in ACS/HRC configuration (dataset id# J9OT07070). 
The Chandra X-ray image was a 0.5-8.0 keV ACIS image from an observation ta ken on 2008 



uly 9 -11. The X-ray image was reprocessed following the procedure described in 



Park et al. 



( 20061) . We removed t he pixel randomization and applied the 'subpixel resolution' technique 



( iTsunemi et al. 



20011 ) to further improve the effective spati al resolution of the image. It 



1972 



Lucy 



1974h 



was then deconvolved with the Lucy-Richardson algorithm (iRichardsonI 
using a monochromatic PSF of 1.5 keV (2.5 keV for the hard band) generated by the 
Chandra Ray Tracer and MARX. The results were smoothed by a Gaussian of O'.'l FWHM. 
Shown in Figure [8] is a comparison between the HST, 36.2 GHz ATCA, 8.6 GHz ATCA, 
and Chandra X-ray images. The coordinates of the HST and Chandra images have been 
shifted such that the optical and X-ray rings are centered on the 2008 Oct 36.2 GHz 
radio ring with an uncertainty of roughly 50 mas in each coordinate; it also implies the 



central optical nebulosity aligns with the fiat spectrum radio component. However this step 
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precludes us from making further judgements about an asymmetry in the radio expansion 
with respect to the optical ring. One thing that is interesting about this plot, however, 
is that the radio images are more circular than the optical and X-ray. If we apply the 
ridge-following algorithm to the flux-unweighted optical ring in the HST image, the width 
of the fitted ellipse semi-major axis is (y.'78 and the semi-minor axis is C/.'S? with about 40 
mas uncertainty in each axis. Compared to the radio ring of size 0'.'79 x 0'.'68 with similar 
uncertainty, the optical ring is slightly more elliptical. The difference in ellipticity can be 
explained by the difference in location at which the different types of emission are generated. 
Non thermal radio emission is produced as the result of stochastic acceleration of particles 
at a shock front, and is primarily dependent on density for particle injection into the shock. 
Thermal X-rays are produced by material heated by the shock and is dependent on both 
density and temperature of the heated material. Hence we would expect radio emission 
to trace shock fronts to higher latitudes above the equatorial plane and X-ray and optical 
images to trace the higher density of heated material in the equatorial ring. 
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Fig. 6. — Results from the ellipse fits to the morphology of the combined 2008 October 
and 2008 April 36.2 GHz observations. The green ellipse drawn atop the image is the flux- 
weighted ellipse fit to the red crosses returned by the ridge-following algorithm. The centroid 
of the green ellipse is shown as a white cross almost hidden under the red diamond which 
represents the centroid of the flux - unweighted ellipse fi t. The yellow star is the measured 



optical position of the progenitor (IReynolds et al. 



19951 ): this moves to the position of the 



small yellow circle when the proper motion of the remnant and the LMC is accounted for. 
The black circle represents the position of the fiat spectral index feature. 
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A Ra (arcsec) Distance from predicted position (arcsec) 

Fig. 7. — Radial profiles overlaid on the combined 2008 Apr and 2008 October 36.2 GHz 
image. Radial lines are propagated away from the predicted progenitor position in (a) and 
their radial distances are measured in (b). The positive distances measured in (b) correspond 
to rays propagating through the left half of (a). 
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In figure [9] the HST and Chandra images have been overlaid with the radio contours 
from the 36.2 GHz radio image. From the peaks in the radio contours it is evident that the 
forward shock is only shghtly larger than the O'.'S optical ring. The truncated-shell models 
fitted to the 2008 October data give a more accurate shell radius of 0'.'85 ± 0'.'06 for the 
36 GHz data which is larger but still consistent. However the radius measured at 8.6 GHz 
is 0'.'892 ± 0'.'002, which suggests that the forward shock has almost completed its transit of 
the optical ring. 
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Fig. 8. — SN1987A at different wavelengths, (a) is the 2006 Dec 6 HST image, (b) is the 
2008 Oct 7-12 ATCA 36.2 GHz continuum image at robust=0.5 weighting, (c) is the 2008 
Oct 9-10 ATCA 8.6 GHz super-resolved image at uniform weighting and (d) is the 2008 Jan 
9-11 Chandra 0.5-8.0 keV image. The 8.6 GHz and 36.2 GHz images have a greyscale range 
of 0-12 mJy beam~^ and 0-4 mJy beam^^. Images (a) and (d) have had their coordinates 
shifted to coincide with the radio images as in Figure [9l The linear greyscale on each image 
has been stretched to highlight intrinsic low brightness features. 
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4. Conclusions 

In this paper we have presented results from the 36.2 GHz extension of the monitoring 
observations of SN1987A performed in April and October 2008. The October observation 
was conducted under conditions of good atmospheric stability. A flux density of 27 ± 6 mJy 
was obtained which allows us to derive a spectral index of a ~ —0.83 when fltted along 



with flux densities from the 
technique of 



Ng et al. 



ower frequency monitoring observations. The Fourier modeling 
(j2008l ) was employed to produce a model flt of the torus that is 
consistent with the 8.6 GHz image. A comparison of the 36.2 GHz image with the 8.6 GHz 
image at similar resolution was used to obtain a spectral index image that showed a region 
of flatter spectral index about 0'.'2 east of the progenitor. Simulations suggest a component 
of this may be an artifact of the imaging process. However follow-up observations are 
suggested to determine if the remaining component is real central region emission not 
fltted by our Fourier models. Higher sensitivity observations with the newly-commissioned 
Compact Array Broadband Backend (CABB) will help further clarify this feature. An 
elliptical flt to the 36.2 GHz radio morphology shows that the centroid of the remnant is 
within positional accuracy of the progenitor's location when corrected for the proper motion 
of the LMC and the motion of SN1987A within the LMC. The flt is also located within 
positional uncertainty of the flat spectral index feature. Comparisons between radio and 
X-ray images show that the remnant is more spherical in the radio, which can be explained 
by two different emission mechanisms at work. From the Fourier models and comparisons 
to the optical emission we conclude that the forward shock in the remnant is completing its 
transit of the optical ring. 
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John Kirk for their insightful comments on particle acceleration at shock fronts. This 
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for non-HST data is provided by the NASA Office of Space Science via grant NAGS- 7584 
and by other grants and contracts. C.-Y.N. and B.M.G. acknowledge the support of the 
Austrahan Research Council through grant FF0561298. 



- 29 - 




Fig. 9. — Overlay of the combined HST 2006 Dec 6 optical (green), Chandra 2008 Jan 9-11 
X-ray (blue), and ATCA 2008 Oct 36.2 GHz radio images (orange-yellow contours) formed 
by shifting the optical and X-ray coordinate systems to center on the radio ring from the 
2008 Oct 36.2 GHz radio image at robust=0.5 weighting. Radio contours are at 14 (orange), 
30, 40, 60, 70 and 85% (yellow) of the maximum at 2.4 mJy beam~^. The outermost contour 
and the contour within the optical ring are at the same 14% level. 
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